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In this paper we study the dynamical behaviour of a simple cosmological model defined by 
ON . a spatially flat Robertson- Walker geometry, conformally coupled with a massive scalar field. We 

determine a Lyapunov-like function for the non-linear evolution equations. From this function we 
prove that all the stationary solutions are unstable. We also show that all initial conditions, different 
from the stationary points, originate an expanding universe in the asymptotic regime, with a scale 
»' | | parameter a(t) that goes to infinity and the scalar field 4>(t) that goes to zero in an oscillatory way 

6J[). . We also find two asymptotic solutions, valid for sufficiently large values of time. These solutions 

correspond to a radiation dominated phase and to a matter dominated phase, respectively. 



H ' I. INTRODUCTION 

In this work we briefly analyze some new properties of an oversimplified cosmological model: a spatially flat 
Robertson- Walker geometry, conformally coupled with a massive scalar field, which has been studied in many papers 
0. The universe is well described by spatially flat (k = 0) cosmological models, which are also compulsory if we use 
inflationary theories, and an inflation field is usually included in these models. Moreover, the conformal coupling is 
necessary if we want to satisfy the equivalence principle ||. Hence, the studied oversimplified structure is necessarily 
contained in the universe, albeit complemented by much more detailed features. Since this simplified stru 

cture by itself implies several of the most important cosmological properties we believe that the results of this letter 
are of interest. 

These properties are: 

a. - It expands to infinity for aH-iion-stationary initial conditions. 

b. - It contains a Lyapunov-like Elfunction F with a positive time derivative and therefore is essentially time asym- 
metric. 



J A Lyapunov function is an evergrowing function of t that vanishes at a critical point and it is different from zero in a 
neighborhood of this point. Our function is not a proper Lyapunov function because it is not of definite sign. 
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c- The final phase of the universe evolution is either radiation or matter dominated, 
d.- The scalar field <p(t) goes to zero in an oscillatory way in the asymptotic regime. 



II. THE MODEL 

Let us consider a flat Robertson- Walker geometry with scale factor a(t) conformally coupled with a neutral massive 
scalar field (j)(t) of mass m. The hamiltonian reads ||: 

where 

k=\{x-y 2 ) + \x 2 y 2 (2) 

x = m(f>, y = ma, p x = x, p y =y and the dots indicate derivative with respect to the conformal time t, defined by the 
relation dt p — a(t) dt , where t p is the commoving or proper time. The field equations read: 

x=-xy 2 , y=yx 2 (3) 

and the Einstein constraint is Tt — 0. This condition expresses the Hamiltonian constraint of general relativity. The 
fixed points are located at: 

x arbitrary y — p x = p y = 

y arbitrary x = p x = p v = 
These fixed points are stationary solutions of the system. 

III. MAIN DYNAMICAL PROPERTIES OF THE MODEL 

Let us define the function: 

F(y, V) = yV (4) 
Taking the time derivative of this function along an orbit of the system and using the second equation (|^) we obtain: 

f =/+^V^0 (5) 

Hence, function F is a monotonous growing function of t. Then, we have 

lim F(y(t), y (*)) = +oo (6) 

t — >+oo 

for an arbitrary non-stationary initial condition. The case lim t ^ +oc F{y(t), y (t)) — C, where C is a constant, is 
excluded since for great values of t we have y(t) = ±\J2Ct + Co and the second equation (||) can not be satisfied for real 
values of x. Here, Co is an integration constant that appears when we solve the differential equation y(t) y (t)) = C. 
We conclude that limt^+oo y(t) = ±oo, i.e. all the orbits in phase space (x,y,p x ,p y ) go to infinity. It is clear that 
eq. (||) forbids the existence of periodic solutions. If a periodic solution x(t) , y(t) of period T exists, integrating both 
sides of eq. (^|) along this solution we obtain zero in the left-hand side and a positive number in the right-hand side 
and we have a contradiction. All initial conditions different from the fixed points yield an universe evolution that end 
in and expanding phase. 

Moreover, from this fact we see that all the fixed points are unstable. We also conclude that the system has no 
chaotic behaviour because all the orbits different of the fixed points are not bounded (see Q ) . 

Let us study now the asymptotic behaviour of the scalar field. For great values of t the function y(t) goes to ±oo 

in a monotonous way, and the product y y is positive. Then, the function G(y, y) = ^ is well defined and positive in 
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the asymptotic regime. Taking the time derivative of G along an orbit of the system and using the second eq. (||) we 
obtain 



The constraint H = gives 



r _ y _y__ V 2 x 2 - V 
y y z y z 



■ 2 .2 9 9 

y = x + x 2 y 2 



(7) 



(8) 



Replacing this expression of y in eq. (W) we obtain 



. 2 

X 

G= — J ^ 



Since the function G is positive and always decreasing we have 



lim G = lim — = k 

t — >+oo t — >+oo y 



(9) 



(10) 



where A; is a non-negative constant. Since function G tends to a constant value in a monotonous way, its first 
derivative must tends to zero when t — > +oo , i.e. 



lim G = lim = 



t — >+oo 



t^ + oo y A 



Equation (H) can be written as follows 



■2 .2 

V _ x 

y 2 y 2 



Taking the limit when t — > +oo in both sides of this equality and using eqs. (llOl) and (111), we obtain 



(11) 



(12) 



(13) 



We will now show that the constant k is zero. From eq. (|13|) we have that limt^ +oc x — ±fc. Let us consider a first 
case where the limit is positive. For sufficiently great values of t, x must approach the value +k in a monotonous way. 
In fact, small oscillatory behaviour with decreasing amplitude around the value +k is impossible, since such type of 
behaviour implies a sequence of changes of sign of the second derivative of x. From the first equation (^) we see that 
it is impossible because the function x does not change sign. 

But, if the function x approaches the value +k in a monotonous way, its first and second derivatives must go to 
zero when t — > +oo. Since lim t ^ +oc y 2 = +oo, we see from the first eq. @ that x can not go to zero. The analysis 
of the case where x approaches the value —k is analogous. Then, we conclude that the constant k must be zero, i.e. 



lim x = k 

t — >+oc 



lim x 







(14) 



We have proved that x can not approach the value ±fc in a monotonous way. The analysis is independent of the 
value of k. Then, we conclude that x has an oscillatory behaviour with a decreasing amplitude when t — > +oo. This 
oscillatory behaviour is not possible for k ^ but it becomes possible in the case k = because x{t) can change sign. 

To summarise, we have proved that the scalar field <fi = 2- goes to zero in a oscillatory way, the factor scale a 2 goes 

to infinity in a monotonous way, and ^ goes to zero when t — ► +oo (cq.[l0] with k = 0). 
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IV. ASYMPTOTIC SOLUTIONS 



For t >> 1 we can obtain explicit asymptotic solutions. Instead of system (y) we will consider the system: 

x=-xy 2 , H = (15) 

because deriving the second equation with respect to t and using the first one we obtain y — yx 2 . System ( |l5| ) 
automatically satisfies the constraint TL — 0, while the system (J3[) has solutions that do not satisfy this constraint. 

Inspired by exhaustive numerical calculations we have been able to determine two different asymptotic solutions 
for t >> 1, namely: 

x(t) ~ j sin Q* 2 J , y(t)~t (16) 

and 

~ j sin Qi 3 ) , y(t) ~ t 2 (17) 

Introducing these expressions in system (O) and taking t >> 1, it can easily be verified that ( |l6| ) and p7| ) are 
asymptotic solutions. Expressions ( |l6|) and (|17|) agree with the general results of section III. 
Expressing these asymptotic solutions in terms of commoving time t p we conclude that: 



a. - The final phase of solution (|16|) is radiation dominated (since yitp) ~ tp). 

b. - The final phase of solution ( |l7j) is matter dominated (since y{tp) ~ ij,). 

We have proved that (|l^) and (|17|) are asymptotic solutions but this does not imply that they are the unique 
solutions for t ^> 1. For a large number of initial conditions chosen at random, numerical computations show that the 
associated asymptotic solution is given by (§!]), i.e. a matter dominated evolution. 

We have performed our analysis in the conformal time t. In general, to translate the results to the commoving time 
t p is not an easy matter (in some cases it is practically impossible) when the scale factor a(t) is not monotonous. In our 
case, for sufficiently large values of t, a(t) is monotonous and the asymptotic behaviour in the variable t can be easily 
translated to the commoving time t p . If we choose physical initial conditions in such a way that a(t = 0) a (t = 0) > 0, 
it is easy to show that a(t) is a monotonous function of t. Then, for such initial conditions all the results obtained 
in terms of t can be easily translated in terms of commoving time t p , for arbitrary values of t. For more general 
cosmological models where the scale factor does not has a monotonous behaviour, the dynamical analysis in terms of 
conformal time t is not justified. In fact, it is clearly wrong from the point of view of dynamical system theory. 

V. CONCLUSIONS 

In spite of the simplified nature of our cosmological model, its dynamical evolution is controlled by highly non-linear 
equations. Nevertheless, we have been able to obtain, in a rigorous way, the most relevant dynamical properties of the 
model: the instability of all the stationary solutions, the expansive nature of the evolution, the fact that the scalar 
field 4>{t) goes to zero in oscillatory way and explicit expressions for the asymptotic solutions, giving the two possible 
phases of the universe. 

Under the time inversion t — > —t the hamiltonian (|^), the field equations (||) and the constraint TL = remain 
invariant, so the system is trivially time-symmetric. Nevertheless, if the motion begins at a finite point with a finite 
value of F (and this would be the case for the real universe that began very small and with high temperature in 
a quantum phase), the motion always goes to infinity. Then, it is time asymmetric since its initial and final states 
are different. Of course, the inverse motion is also a solution of the evolution equations. But the motion towards 
infinity is usually expanding and with creation of matter-field energy, and therefore corresponds to what we see in 
the observable universe. 
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